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Abstract 

We have examined the local structure of PbMgi/gNba/sOs-PbTiOs (PMN-PT) and 
PbZn]^/3Nb2/303-PbTi03 (PZN-PT) solid solutions using density functional theory. We find that 
the directions and magnitudes of cation displacement can be explained by an interplay of cation- 
oxygen bonding, electrostatic dipole-dipole interactions and short-range direct and through oxygen 
Pb-B-cation repulsive interactions. We find that the Zn ions off-center in the PZN-PT system, 
which also enables larger Pb and Nb/Ti displacements. The off-centering behavior of Zn lessens 
Pb-B-cation repulsion, leading to a relaxor to ferroelectric and a rhombohedral to tetragonal phase 
transition at low PbTi03 content in the PZN-PT system. We also show that a simple quadratic 
relationship exists between Pb and B-cation displacements and the temperature maximum of di- 
electric constant, thus linking the enhanced displacements in PZN-PT systems with the higher 
transition temperatures. 
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Heterovalent ferroelectric perovskite solid solutions have been intensively studied [? ? ? 
? ? ] as they exhibit a range of interesting structural and dielectric properties which make 
them useful in device applications. Relating the macroscopic properties of these materials 
to the microscopic properties of the constituent atoms is an important scientific goal and is 
vital for the rational design of new materials with improved properties. The technologically 
important [? ] PbMgi/gNbs/sPs-PbTiOg (PMN-PT) and PbZni/gNba/aPa-PbTiOa (PZN- 
PT) solid solutions are good subjects for an investigation of structure-property correlations 
due to the availability of a wide range of experimental data and because the small difference 
between Mg and Zn ionic sizes [? ] cannot explain the differences in the characteristics of 
the two solid solutions such as a higher and dielectric constant observed for PZN [? ] . 

Ab-initio Density Functional Theory (DFT) calculations were performed for the (1- 
x)PMN-xPT and (l-a;)PZN-2;PT systems using 2 x 2 x 2 40-atom or 3 x 2 x 2 60-atom 
supercells with periodic boundary conditions at experimental volume [? ]. The energy of the 
system was evaluated using a local density approximation exchange-correlation functional [? 
] and was minimized with respect to the atomic coordinates by a quasi-Newton method with 
no symmetry imposed. A 2 x 2 x 2 fc-point sampling of the Brillouin zone was used. The 
calculations were done with designed non-local norm conserving pseudopotentials [? ][1] We 
study four compositions with x = 0, 0.25, 0.625 and 1.0 For pure PMN and PZN we used 
a 60-atom supercell with B-cations arranged according to the random-site model [? ]. A 
B-cation arrangement consistent with the random-site model as well as two disordered B- 
cation arrangements were used to model the 0.25 composition and three disordered B-cation 
configurations were used for the 0.625 composition. 

The excellent agreement between the pair distribution functions (PDF) obtained from 
our relaxed 60-atom structure for PMN and by neutron-scattering experiments [? ] (Fig- 
ure 1) and similar agreement between DFT and experimental PDFs in previous stud- 
ies of PbScs/aWi/sOs-PbTiOg (PSW-PT), PbScs/sWi/gOg-PbZrOg (PSW-PZ) [? ] and 
PbZri_j;Tia;03 (PZT) [? ? ] solid solutions indicate that the size of our supercells is 
sufficient for capturing the local structure of these materials. 

We find that at all compositions PZN-PT and PMN-PT local structure is governed by 
three main interactions. First, bonding interactions with oxygen favor large Pb and B-cation 
displacements (Table I) and creation of short Pb-0 and B-0 bonds. Second, short-range 
repulsive Pb-B-cation interactions create anisotropy in the potential energy surfaces felt by 
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the Pb cations, leading to a variation in the preferred local displacement direction for the Pb 
atoms (Table I, Figure 2). Third, electrostatic interactions favor dipole alignment and tend 
to minimize the scatter in displacement direction, eliminating displacements opposite to the 
overall polarization of the supercell (Figure 2). Breathing motion by the B-Og octahcdra (to 
achieve preferred octahedral volume) and small octahedral rotations (0-5°) are also present. 

Examination of the data in Table I shows that for both solid solutions Pb displacement 
magnitudes predicted by DFT are unchanged by the addition of PT, whereas the average 
magnitude of the B-cation displacements steadily increases with PT content. The cation 
displacements are larger in PZN-PT than in PMN-PT, with especially large difference for 
the B^+ cations. While Mg ions are only slightly displaced, the Zn ions exhibit significant 
distortions even in pure PZN and Ti-like 0.27 A distortions in Ti-rich PZN-PT. As will be 
explained later, the large Zn distortions also account for the larger magnitude of Pb, Nb 
and Ti distortions in the PZN-PT system. 

The contrast between the absence of Mg off-centering and the large Zn displacements is 
due to the different electronic structure of the two ions. While Mg is a simple metal and 
]V[g^+-0 bonding is essentially ionic, the imperfect screening of the nuclear charge by the 
(i-electrons and the presence of low-lying p orbitals make Zn more polarizable and enable 
covalent bonding with oxygen atoms [? ? ]. 

As in other Pb containing ferroelectric perovskite solutions, Pb displacements are the 
determining factor for the average structure of our solid solutions [???]. Inspection of 
Figure 2 shows that the location of Mg atoms is the primary influence on the direction of Pb 
distortions. In all cases Pb atoms avoid cube faces with three Nb atoms and one Mg atom 
and move toward cube faces with two Mg and two Nb atoms. The avoidance of Nb-rich 
faces is due to the presence of oxygen atoms with two Nb neighbors. Such oxygen atoms 
have higher B-0 bond order than the oxygen atoms in PbTiOs due to the higher valence 
of Nb. Since the total bond order of cation-oxygen bonds is conserved at 2, this prohibits 
the formation of short Pb-0 bonds of high bond order and creates a strong effective Pb-Nb 
repulsive interaction. On the other hand, the oxygen atoms with Mg and Nb neighbors have 
a lower total B-0 bond order than the oxygen atoms in PbTiOs, necessitating a formation 
of more short Pb-0 bonds to compensate for B-0 bond depletion [? ? ] and leading to a 
weaker effective Pb-Mg repulsion. 

The larger Pb and B-cation displacements in the PZN-PT system are due to the couphng 
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between Zn and Pb and Nb off-center displacements through the Pb-B-cation repulsive 
interactions. The ionic size of Mg and Zn atoms is essentially the same [? ], which leads 
to very similar shortest allowed Pb-Mg and Pb-Zn distances. The greater ability of Zn 
atoms to off-center allows Pb atoms in PZN-PT solution to preserve the required Pb-B- 
cation distances for larger Pb displacements than in the PMN-PT system. The larger 
Pb displacements in turn cause larger Nb and Ti distortions to prevent high Pb-B-cation 
repulsion at short Pb-Nb and Pb-Ti distances. 

The electrostatic dipole-dipole interactions favor alignment of cation displacements and 
prevent cation displacements in the (100) direction in direct opposition to the overall polar- 
ization along [100] in Figure 2. The local direction of Pb displacement is determined by a 
balance between the need to align with other cation displacements and the need to minimize 
direct and through-oxygen Pb-B-cation repulsion [? ? ]. 

At high Ti content, coUinear [100] Pb atom distortions create short Pb-0 bonds, maxi- 
mize dipole alignment and minimize local Pb-B-cation repulsion by avoiding the Ti cations 
located along the [111] direction. As Ti content decreases, the Pb local environments become 
anisotropic due to the presence of the overbonded and underbonded oxygen atoms which 
block Pb distortions toward some [100] perovskite cube faces. This leads to a transition to a 
rhombohedral phase, where Pb cations move in a variety of low symmetry directions around 
the overall [111] polarization direction toward the available low-repulsion cube faces. With 
the onset of random-site B-cation ordering, the high anisotropy created by a large amount of 
faces blocked by the presence of overbonded oxygen atoms forbids the formation of a long- 
range ferroelectric state even in the [111] direction and leads to a relaxor phase. A similar 
explanation, albeit focused on the B-cation size difference, for the relaxor to ferroelectric 
transition in the ordered PMN-PSN system, was proposed by Farbcr and Davies [? ]. 

While exact Pb-environment population analysis in PMN-PT and PZN-PT is currently 
not possible due to the complex and unknown nature of the short-range B-cation ordering 
in these systems [? ? ], the above framework coupled with the differences in Zn and Mg 
distortion behavior explains why in the PZN-PT system the tetragonal phase is present at 
lower PT content (10%) than in the PMN-PT system (35%), as well as the recent discovery 
of ferroelectric phases even in pure PZN itself [? ]. A large Zn distortion softens Pb-Zn 
repulsion, making perovskite faces with Zn atoms more friendly toward a Pb distortion than 
the corresponding perovskite faces with Mg atoms in PMN. This increases the population of 
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the local environments favorable to tetragonal distortion, moving the relaxor to ferroelectric 
transition and the morphotropic phase boundary (MPB) to lower Ti concentrations. 

We find that there is also a strong relationship between the low-temperature cation dis- 
placements as found by our DFT calculations and the experimentally obtained temperature 
position of the dielectric constant maximum T^^max- A qualitative relationship between av- 
erage lattice displacement and ferroelectric to paraelectric phase transition temperature Tc, 
which is closely linked to Te^max, was found by Abrahams, Kurtz and Jamieson (AKJ) in 
1968 [? ]. However, examination of the original AKJ paper shows that the largest dis- 
crepancies between predicted and experimental Tc are for the Pb containing perovskites 
PbFei/aNbi/aOg and PT. We find that T,^^^^ at 1 MHz for PMN-PT, PZN-PT, PZT, 
PbSci/sNbi/aOg (PSN), PSW-PT and PSW-PZ solid solutions is predicted by 



Te,raax — ^C^Pb + ^C^B'/b'; (l) 

where cipb is the average magnitude of the Pb distortions, d^' is the average magnitude of 
the distortions, /b' is the fraction of the ferroelectrically active B-cations in solution and a 
and b are constants (1739 and 5961 respectively) in units of K/A^. 

The dependence of T^^max on local structure presented in Equation 1 can be further 
simplified by transforming cation and oxygen displacements Uj into K local dipole moments 
Poi via their Born effective charge tensors Z* This leads to 




where 7 is a constant, is the number of primitive unit cells and Vq is the volume of a 
primitive unit cell and |Poi| and Pq cire the K average magnitude of the local dipole moment 
and the overall polarization respectively. Figure 3 shows the correlation between the max 
predicted by Eqns. 1-3 and experimental T^^max data [? ]. The clustering of the data around 
the y=x line indicates that Eqns. 1-3 accurately and quantitatively capture the trends in 
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Te,max among these fairly different systems. Comparing tfie data for PMN-PT and PZN-PT 
systems in Table I, we see that the higher Tg^max for PZN-PT system is directly related to 
larger Pq and |Poi| values due to off-centering behavior on the site. 

The correlations presented in Figure 3 can be explained based on fourth-order Landau 
theory for ferroelectric crystals. The equilibrium polarization is given by 

dG/dP = = 2a{T-To)P + 4(3P^ (4) 

where a and /S are constants, Tq is the Curie-Weiss temperature at which the dielectric 
constant e diverges and P is polarization. Substituting T=0 used in our DFT calculations 
and solving for Tq, we obtain Tq — 2/9Pq/q;, thus relating the Tq to the square of computed 
K polarization. The high quality of the fit in Figure 2 means that for the Pb-based 
perovskites examined here the ratio (5 /a is nearly constant. 

The energy difference between ferroelectric and paraelectric phases (the ferroelectric in- 
stability AE'fe) is commonly thought to correlate with Tq. However, AE'fe does not correlate 
with To for Pb-based perovskites such as PZT [? ], PSW-PT [? ] and (Cd,Pb)Ti03 [? ? ] 
(CPT). The constancy (3/a helps explain this observation. In fourth-order Landau theory 
A£;fe at T = is given by A£;fe = oi^T^/Ap, with T^^{AI3/a^)l^EY^. The (3 parameter 
scales the leading order repulsive term and corresponds to the repulsive potential felt by 
the off-centering Pb and B-cations. This potential varies from system to system based on 
different sizes of the B-cations [? ? ]• Since the ratio of j3/a is constant for a wide variety 
of ferroelectric solutions as shown above, the value of varies. This makes A£^fe a poor 
predictor of Tq. 

The nearly constant dependence of Tg^max on Pq and |Poi| can be interpreted as follows. 
Experimental results have shown that ferroelectric to paraelectric phase transitions in Pb 
based perovskites such as PMN and PZT are of order-disorder character [? ], with Pb 
displacements at Tq of up to 70% of the low temperature value. The parameter governing 
the temperature of such an order-disorder transition is the coupling strength between the 
Pb and B-cation displacements in neighboring unit cells (corresponding to the J parameter 
in spin models). The coupling is due to electrostatic effects favoring dipole alignment and 
covalent chemical bonding effects disfavoring under- and over-bonded oxygen atoms. Dipole 
alignment strengths depend only on polarization strengths and gives rise to the overall 
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dependence of To on Pq or |Poi| , with the variation in through-oxygen couphng accounting 
for the differences between experimental and predicted T^^max values in Figure 3. While 
predicted T^^max values arc of the same quality for Eqn. 1 and Eqn. 2, the correlation between 
Pq and Te^max is somewhat worse. This suggests that the freezing in of large distortions 
enhances stability of ferroelectricity, giving rise to higher Tg max- Even though some of the 
distortions frozen in are not parallel, they still contribute to the coupling of dipoles in 
neighboring unit cells due to through-oxygen interactions. 

In summary, we have examined structure-property correlations for PMN-PT and PZN- 
PT solid solutions. We find that the local direction of Pb displacements in these systems 
is governed by an interplay between the electrostatic dipole-dipole alignment interactions 
and Pb-B repulsive interactions, as was previously found for PZT. We find that Zn ions 
off-center; this reduces Pb-B repulsions resulting in smaller Ti content at MPB and allowing 
larger displacements by the Pb, Nb and Ti ions. A strong relationship exists between 
the cation displacements or Pq and Te max, which can be explained in the Landau theory 
framework by the constant ratio of a and /? parameters and is a natural consequence of the 
order-disorder nature of phase transitions in Pb-bascd fcrroelectrics. The link between a and 
(3 also explains the opposite trends in A£;fe and found in the PZT, PSW-PT and CPT 
systems. Fundamentally, we have shown that behavior of individual ions in heterovalent 
perovskite ferroelectric solid solutions can be directly linked to changes in local structure 
and macroscopic collective properties. 
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TABLE I: Results of DFT calculations for PMN-PT and PZN-PT systems. Pb displacement 
magnitudes in PMN predicted by DFT are in agreement with recent experimental data [? ] and 
P magnitudes for PMN and PT are in good agreement with previous theoretical calculations [? 
]. Predicted cation displacements from center of oxygen cage in A, Pb displacement angle scatter 
^Pb in °, average and local polarization in C/m^ and experimental T^^max [? ] in K. 

Pb B2+ Nb/Ti^pb Po T„ 
disp disp disp avg,loc exp. 

PMN 0.398 0.056 0.168 65 0.38, 0.67 276 

PMN-0.25PT 0.389 0.080 0.181 33 0.55,0.68 397 
PMN-0.625PT 0.387 0.099 0.216 28 0.66, 0.74 583 
PT 0.440 0.280 0.88,0.88 765 



PZN 0.444 0.148 

PZN-0.25PT 0.461 0.258 

PZN-0.625PT 0.424 0.270 

PT 0.440 



0.189 67 0.43,0.73 424 

0.205 35 0.66,0.79 547 

0.237 27 0.74,0.80 643 

0.280 0.88,0.88 765 
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FIG. 1: Comparison of PDFs obtained from relaxed structure of 60-atom supercell of PMN (solid) 
and by neutron-scattering (dashed) [? ]. Similar agreement is obtained for other systems [? ? ]• 



[1] Ionized +2 reference configurations were used for Pb, Mg, and Zn, ionized +4 and +5 were 
used for Ti and Nb and a neutral reference configuration was used for O. Energy cutoff of 50 Ry 
were used for all elements except for Zn for which 56 Ry was used. Real space cutoffs were 2.0, 
1.83, 1.95, 1.70, 1.75 and 1.5 for Pb, Mg, Zn, Ti, Nb and O respectively. 

[2] Following approach of Ref. 26, we used approximate Z* values of 3.6, 2.6, 2.6, 7.4, 4.0, 12.0, 
6.4 and 6.4 for Pb, Mg, Zn, Nb, Sc, W, Ti and Zr respectively. Oxygen Z* were ~ -5.0 and 
~ -2.5 for the parallel and perpendicular elements of the tensor respectively, varied slightly to 
satisfy the electrostatic sum rule for each composition. 
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FIG. 3: Correlations between experimental T. ^ax and predicted T,,max for PZN, PZN-PT, PMN, 
PMN-PT, PSW, PSW-PT, PSW-PZ, PSN and PZT systems. Values for predicted T.^max are 
obtained by Eqn. 1 (squares) with a=1739 and b= 5961, Eqn. 2 using local polarization with 
7=942 (stars) and Eqn. 2 using overall polarization with 7=1189 (circles). 
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